Magnetic quantum tunneling in the half-integer high spin 
molecule CrNie : a ^SR study 
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In high spin molecules metal ions are coupled by ferro or antiferromagnetic short range 
interactions so that their magnetic moments are parallel or antiparallel to each other at 
temperatures (T) much smaller than the coupling constant J. This leads to a high spin 
(S) value in the ground state. In the crystal lattice, the molecules are well separated from 
each other and the active part of the Hamiltonian at ksT <C J is TC = —DS 2 , so that 
up and down spins (S z = ±S) states have identical energies, and S z can escape from one 
state to the other via either over-the-barrier motion, or tunneling. So far, the escape rate 
T in HSM was found to have smooth T dependence throughout the cooling process from 
DS 2 <C kgT to k^T <C DS 2 P]J2|]. In this case, the transition from the over-the-barrier 
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motion to tunneling, is referred to as second-order transition ||. In the present work, 
we use the muon spin relaxation technique to show evidence for temperature-independent 
tunneling at ksT <C DS 2 , in zero applied magnetic field (ZF), in the HSM CrNig. Moreover, 
the crossover from the classical to quantum regime is sharp, and appears to be consistent 
with a first order transition of the escape rate |3j . 

In a first order transition, over-the-barrier motion is abruptly replaced by tunneling 
between ground states, as the temperature is lowered down to the barrier height. 

The CrNi6 molecule is described in details in the methods section. The magnetization 
M vs magnetic field H measured at 6 K in this compound is depicted in Fig. |l|. M displays 
a saturation value of 15 //b, arid fits well to a Brillouin function of S = 15/2. This S value 
is expected from a Cr(III) (Sc T = 3/2) ion ferromagnetically coupled to 6 Ni(II) (Sn; = 1) 
ions. The ferromagnetic nature of the interaction is easily rationalized in the frame of a 
localized electron model, since the quasi-linear Cr-CN-Ni unit allows orthogonality between 
the Cr(III) and Ni(II) wave functions. 

Also shown in Fig. |I] is the thermal variation of the molar susceptibility as the product 
XmT vs T, which monotonously increases from room temperature to 6K. The fit of the 
measured susceptibility to the one expected from the Hamiltonian H = — JJ2m ^Cr • Sni — 
D(Sl) 2 , where S* is the total spin of the molecule, gives a good agreement when choosing 
a mean g value = 2, a J value « 24K, and \D\ (S 1 ) 2 = 4.2 K | (see Fig. |l|). The value of 
| D | S 2 is in agreement with ac-susceptibility measurements ||. S 1 ' = 15/2 is well separated 
in energy by 3 J/2 (36K) from an excited level where S l = 13/2 , itself 3 J/2 below a second 
level where S l = 11/2, etc. 

The /iSR experiments presented here were performed in PSI after preliminary investiga- 
tion at ISIS using the longitudinal field (LF) technique. In this technique one follows the 
time evolution of the spin polarization P z (t) of a muon implanted in a sample, through the 
asymmetry A(t) oc P z if) in the positron emission of the muon decay 0. In addition, an 
external field is applied along the initial muon spin (longitudinal) direction which defines 
the z axis. In view of our experimental results, it is worth discussing the expected behavior 
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of P z {t) in two cases: (I) the field at the muon site is random and static or (II) the field is 
dynamically fluctuating. In the static case, and when = 0, the muon polarization relaxes, 
and the time scale for relaxation [r] is determined by ( / y fl {B])~ 1 where 7 M = 85.162 MHz/kG 
is the muon gyromagnetic ratio and [B] is the scale of the local magnetic field. However, 
not all muons will relax since there is a statistical probability that some of them (effectively 
1/3) will reside in a site where the field is parallel to the muon spin. Therefore, after a 
minimum in P z (t), we expect lim^oo P z (t) = 1/3 , a phenomenon called recovery |J. When 
the external field is applied, and [B] <C H^, the vector sum of and the internal field 
is nearly parallel to the muon spin, and we expect lim f ^ 00 P 2 (t) = 1. This phenomenaon 
is called decoupling 0. In the dynamical case the muon relaxation is determined by the 
magnetic spectral density at the frequency u>l = l^H^ corresponding to the muon Zeeman 
energy levels (see below). 

In Fig. ^|we show the asymmetry, obtained in CrNig, in ZF, T = 50 mK, and at various 
values of H^. In the ZF case we find lim^oo A(t) = 0. This limit suggests that the field at 
the muon site is not static. Moreover, the time scale for the muon relaxation is 1 //sec. If 
this relaxation was due to static fields it would imply a local internal field in the order of 
tens of Gauss. Such a small field would have decoupled [A(t) = A(0)} with H L = 2 kG or 
higher, in contrast to the observed behavior. We thus conclude that the CrNig molecules 
are dynamically fluctuating even at T = 50 mK and in zero applied magnetic field. We 
interpret these fluctuations as tunneling within the ground spin state, between S z levels at 
the bottom of the anisotropy well, since at 50 mK only these levels are populated. A similar 
conclusion was reached by Vernier et al. using high frequency experiments. 

In order to illuminate this conclusion we show in the inset of Fig. |2| a similar experiment 
in a powder of the Ising spin glass Fe .o5TiS 2 (T g = 15.5 K) at T = 4 K 0. Here, as in the 
HSM, the most important degree of freedom is S z , and the classical ground state is made 
of spins randomly oriented in either up or down directions. However, in contrast to CrNig, 
we can clearly see in Feo.o5TiS2 a recovery of the asymmetry. In addition, the time scale of 
muon relaxation in Feo.o5TiS2 is an order of magnitude faster than in CrNig: yet a field of 
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720 G is sufficient to nearly fully decouple the muon polarization. Thus, Feo.osTiS2 fulfills 
all our expectations from a static magnet and emphasizes in contrast the dynamical nature 
of CrNi 6 . 

Next we determine the fluctuation rate v (proportional to the escape rate T), where 
v is defined by (S t (t)S t (0)) = ((S'*) 2 ) exp(— ut). Since at low temperatures the spins on 
the molecule are locked to each other, the fluctuation rate of an individual spin should be 
identical to that of the entire molecule. Therefore, the relaxation of the muon spin, which 
is mostly coupled to the electronic spin closest to its stop site, will be determined by v. 
A second important parameter is the root mean square of the field at the muon site A 
(= 7 M (Bf) 1 ^ 2 ) where i is a spatial direction. We notice that the muon spin relaxation has 
a square-root exponential shape, which could be explained by a distribution of A, due to 
a distribution of muon sites, according to p(A) = ^Jljix-^ exp ^— Usually, one can 
separate the two contributions to the muon relaxation rate by using the LF data. Here the 
situation is more intricate as the impact of the applied field on the tunneling rate is not 
evident and in a powder sample H L also couples to S x and S y , and induces tunneling, or to 
S z , and reduces tunneling (for most field values). In the absence of any model, we assume 
that v is field independent for the fields (gfisSH < DS 2 ) used here, so as to get an estimate 
of v. These assumptions lead to 



where, A Q is the initial asymmetry, and A^, represent muons that missed the sample ||. 



find reasonable agreement between theory and experiment with a fluctuation rate v = 85(5) 
MHz and a = 10.0(5) MHz (~ 100 G). The fit, shown in Fig. ^| by the solid lines, suggests 
that indeed the external field might have only minimal impact on the powder averaged 
tunneling rate at base temperature. One possible explanation for this is that the interaction 
responsible for tunneling is much stronger than the Zeeman coupling to the external field. 
Now we discuss the temperature dependence of the muon relaxation. In Fig. |3] we depict 




(1) 



We fit our data to Eq. p] using A , Ab, and u, as global parameters for all data sets. We 
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the muon asymmetry in LF of 100 G and various temperatures; the small field of 100 G 
is applied in order to decouple the contribution from nuclear moments to the relaxation. 
We see that from T = 50 to 6 K the muon relaxation rate increases, but between 6 K and 
50 niK it remains constant (no change in the asymmetry). The fact that the transition 
temperature is ~ 6 K reinforces our basic assumption that the electronic spins determine 
the muon relaxation, since this temperature is of the same order of magnitude as the energy 
barrier. The T independent fluctuation rate below ~ 6 K implies that these fluctuations 
stem from QTM, as opposed to thermally activated, over-the-barrier motion. 

We fit the data obtained at different temperatures to Eq. [I], with v the only free parame- 
ter, and plot the result in Fig. |. Clearly, the temperature dependence of v changes abruptly 
at T ~ 10 K. Also shown in this figure by a solid line is a fit of v to an Arrhenius-like law 

u = u q + u c exp(-U/T). (2) 

where v q = 85 /isec -1 . We find v c = 12.0(3) nsec -1 , and U = 52(3) K. The large value 
of U indicates that for T > 10 K the total spin states S f = 13/2 and S f = 11/2, become 
populated (recall that Jm-Cr = 24 K and that the gap between the lowest total spin states 
is 36 K). A comparison of Fig. ^ with Fig. [I| suggests that between ~ 10 K and ~ 6 K the 
molecules condense in a total spin state S* = 15/2. The behavior below ~ 6 K suggests that 
once the temperature becomes as low as the barrier height (\D\ {Si) 2 = 4.2 K), tunneling is 
taking place between the S z = ±15/2 states; otherwise, it would have shown temperature 
dependence down to ~ \D\ [(15/2) 2 — (13/2) 2 ] = 1.05 K. This behaviour appears closely 
related to a first order quantum-classical transition of the escape rate ||. It should be 
pointed out that the nature of the transition is inferred from a nearly zero field experiment 
and is independent of the absolute value of u, which was determined using the LF. 

Our result renders CrNi 6 the first example of a high spin molecule where the tunneling is 
temperature independent at low T, and the transition from the classical activated behavior 
to the quantum one is sharp and consistent with first order. 
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Methods 

CrNi 6 -chemical formula [Cr{(CN)Ni(tetren)} 6 ](C10 4 )9- is obtained by reacting solutions 
of 0.1 molar potassium hexacyanochromate(III) and 0.15 molar Ni(II) tetren perchlorate in a 
mixed solvent water/acetonitrile. The solution is left to evaporate. After a few weeks violet 
parallelepipedic single crystals appear which are collected, filtered and dried in air. Single 
crystal X-ray diffraction shows that the compound crystallises in the monoclinic system, 
P21/n space group, cell parameters a = 15.375 A, b = 24.280 , c = 16.141 A, a = 90°, 
(3 = 115.580°, 7 = 90°. Disorder on the perchlorate anions and on the carbon framework of 
the terminal amine ligands impedes a complete determination of the structure. Nevertheless, 
it appears clearly that the system is made of heptanuclear entities where the Cr(III) ion is 
surrounded by six cyanide ions, each bonded to a Ni(II) tetren (see Fig. [1]). The coordination 
sphere of Cr(III) and Ni(II) can be described as a slightly distorted octahedral. The Cr- 
CN-Ni distance is 5.23 A. The Cr-CN angle is 176.02°; the CN-Ni angle is 170.46°. The 
heptanuclear units are far away from each other (nearest Cr-Cr distance = 14.77 A and 
nearest Ni-Ni distance from neighbouring units = 8.84 A). 
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FIGURES 

FIG. 1. Thermal dependence of the molar susceptibility xmT vs T, and Magnetisation vs H at 
6K (symboles). The fits (lines) are described in the text. Also shown is the CrNig molecule. 

FIG. 2. Muon asymmetry verseus time in CrNig at T = 50 mK and various applied longitudinal 
magnetic fields. This figure demonstrates the dynamical nature of the electronic spins (see text). 
For comparison the inset shows the muon asymmetry in the Ising spin glass Feo.o5TiS2 (T g = 15.5 K) 
at T = 4 K where these spins are static. 

FIG. 3. Muon asymmetry versus time in CrNig at = 100 G and various temperatures 
(symbols). The solid line is a fit to Eq. 1. 

FIG. 4. Temperature dependence of the molecular spin fluctuation rate. The solid line is a fit 
to Eq. 2. 
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